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Abstract 

CP asymmetry in some exclusive decays of charged B meson are calculated in 
a Bethe-Salpeter approach. Hadronic final state interactions are ignored. Com- 
plex decay amplitudes are assumed to arise entirely from perturbative quark- 
antiquark loops. Calculations are done both with and without the gluon quark- 
antiquark vacuum polarization loops. The effects of neglecting the imaginary 
parts arising from the diagonal quark-antiquark loops are also studied. 
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1. Introduction 



It has long been recognised that CP asymmetry in partial decay rates can be gen- 
erated in processes that involve (i) at least two different Cabibbo-Kobayashi-Maskawa 
(CKM) angles and (ii) at least two different strong interaction phases|I|]. These strong 
phases could be generated either perturbatively through the penguin processes or 
through hadronic final state interactions (fsi). In processes involving a single isospin 
in the final state, CP asymmetry is still possible provided that the penguin processes 
generate strong phases 0. 

The commonly followed practice in evaluating decay matrix elements has been to 
use the factorization approximation which reduces the computation to the knowledge of 
certain form factors. This method treats g^, the invariant momentum transfer carried 
by the gluon or the photon in the strong or electroweak penguins, as a free parame- 
ter whose value is usually chosen in the range < < Tnl/2 for B decays. In 
B decays, perturbative strong phases generated by cc states in the penguin loops are 
particularly sensitive to the choice of g^. In contrast to the factorization scheme where 
one is forced to choose an effective value of g^, in a complete calculation would be 
integrated over and would not remain a free parameter. 

In the past, Simma et.al.[0 have employed the formalism of Lepage and Brodsky[^ 
to calculate CP asymmetries in inclusive and exclusive decays of charged B meson. In 
their approach, due to the fact that extra gluons have to be used to accelerate the spec- 
tator quark, new absorptive parts are generated. They also found that absolute partial 
rates were too low by as much as two orders of magnitude. They had to scale up the 
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rates so that the semileptonic rate B DW was correctly normahzed. It is our view, 
as shown later, that it is doubtful if a single process-independent scale factor would 
be adequate for all exclusive channels. The advantage of their method, aside from 
pedagogy, is that it goes beyond the factorization approximation and, consequently, 
is not left free parameter. 

In our work, we have used the Bethe-Salpeter (B-S) method of Mandelstam]^] and 
Nishijima|p to treat the bound state mesons. This method has been applied previ- 
ously to other decay processes, including radiative decay |]^. Before embarking on the 
discussion of our calculations, we wish to state that our results on CP asymmetries 
are somewhat different in size, though not in sign, to those of 0. This is not un- 
expected considering the differences in methodology and approximations used. The 
results, as emphasized in[Q, are trustworthy to within a factor of 2. However, we share 
withQ the diffculty in getting the absolute rates right. We also need to normalize the 
calculated rates such that the semileptonic rates are correctly normalized. This, how- 
ever, does not affect CP asymmetry as absolute normalizations cancel in the definition. 

2. Preliminaries 

We work with the following effective Hamiltonian for h ^ s transition (for h ^ d 
transition, replace s by d)|^ 



where Vqi, etc are the CKM matrix elements, Ci are the Wilson coefficients (discussed 
below) and Oi are the following Current ® Current operators: 




(1) 
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Tree level operators: 

01 = {sq)v~A{qh)v-A, Ol = {Saqi3)v-A{qpha)v-A\ (2) 
QCD penguin operators: 

O3 = {sh)v~AY^{qq')v-A, Oa = {Sabf3)v-AEq'{q'(3q'a)v-A, 
q' 

O5 = isb)v-Aj2(^'l')v+A, Oq = {Sabf3)v-Aj2q' iq'pq'Jv+A] (3) 

Electroweak penguin operators: 
3 

O7 = -(s6)y-A$^ (eg'g'g')v+A, Og = |(s„6/3)y_^ Eg' (eq'q'pq'Jv+A, 
q' 

3 

O9 = ^isb)v-Aj2i'^i'^'l')v-A, Oio = ^Sabf3)v-Aj2q' {eq'q'pq'Jv-A (4) 
q' 

where the subscripts V ± A represent 7^(1 ± 75), a and P are color indices. J2q' 
represents a sum over u,d,s and c quarks. 

Wilson coefficients, Ci, are renormalization scheme dependent when the next -to- 
leading order corrections are considered. However, as the physical quantities should 
be scheme-independent, one needs to define effective Wilson coefficients to be used in 
computation. The manner in which these effective coefficients are defined is detailed 
in 1^, ^, 10|. We only summarize the final result: The effective coefficients are 
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where Nc is the number of colors and 

Ci = 1.1502, C2 = -0.3125, C3 = 0.0174, C4 = -0.0373, C5 = 0.0104, Cq = -0.0459, 
Cy = -1.050 X 10-^ Cs = 3.839 x 10"^ Cg = -0.0101, Cw = 1.959 x 10^=^. (6) 
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To one-loop approximation, the contribution to Ps,e from penguin diagrams with 
insertions of tree operators Oi^2 (see Fig. 1(a)), are 



P^'' = ^CM'^ + hn^-G{m„f.,q% (7) 
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aem(/i) r^ , , , Clifi) 10 2 m 
Pe = [C2{fi) + ^;^][— + - log -^-G{mg,fi,q )]. (8) 

Here the flavor q stands for u and c in the penguin loop and is the invariant 
momentum carried by the gluon or the photon. G{mg, /i, g^) is given by 

G(mg, /i, g^) = — 4 / dxx{l — x)\n[l — x{l — x) — -], (9) 

J 

Obviously, as g^ > 4mg, a strong perturbative phase is generated by virtue of G{mg, /i, g^) 
becoming complex. This perturbative phase arising from the complex nature of G{mq, /i, g^) 
has nothing to do with the strong phases generated by long range fsi among hadrons. 
The latter respect isospin symmetry, while gluons operating in color-space, the pertur- 
bative phases have nothing to do with isospin space. 



In order to accommodate some of the cancellations required by the CPT theoremjll 



we also consider a 1-loop correction to the gluon propagator. Fig. 1(b), leading to a two- 
loop expression for Pg, 



^(,) ^ l_^a^^,^^^^^^10^2_ 1^ ^2)] X: In ^ - G{m,,f^, q% (10) 

It is important to note that we will not concern ourselves with the soft hadronic 
fsi; all phases are considered to be generated by the complex value of G{mq, /i, g^). 



3. Calculation of hadronic matrix elements 
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The Feynman rules for the B-S formahsm of Mandelstam|^ and Nishijima^ have 
been abstracted in Esteve et.aL0. We illustrate the application of those rules by 
writing down the matrix element for D^D^ transition. In Fig. 2, we have shown 

the tree and the penguin graphs 

(11) 
(12) 

X^;^-(^-?+Y)7^(1-75)]. (13) 

Here Ci{q^) is a particular Wilson coefficient which may carry a dependence via the 
function G{mg, fi, q^) and could be complex. Notice that difference between Figs. 2(b) 
and 2(c) is that the penguin is "time-like" or "space-like", in quark language, respec- 
tively. In the B-S formalism, both of these graphs can generate complex amplitudes 
as q is an integration variable spanning the entire space. F'^ is the appropriate matrix, 
(V±A), depending on which of the Wilson coefficients is chosen for Cj(g^). Xb{s) etc. 
are the phenomenological B-S relativistic wavefunctions which we choose according to 
Esteve et.al. 0, 

XM(g) = iV75(l + 5/P)$(P,g), x = 7oX+7o (14) 
where the momentum wavefunction $(P, g) is, 

$(P, q) = exp{--i^[2(^)2 - q']}. (15) 

This is a particularly simple form which in the rest frame of M reduces to 
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^^^^ ^^ restframe ^ ^^^^ 

Clearly, this form assigns the same strength to the space- and time- oscillators. There 
is no compelling reason for this. A more general form could be obtained by assigning 
a different a to (P ■ g)^ and in the choice of g), but we don't concern ourselves 
with such finer details. The values of the parameters 5, a and N are fixed according to 
Esteve et.al. 

5 
a 
N 

nia, mh are the current masses of quarks making up the meson M with the decay con- 
stant /m- 

4. Results 

We present our results under the following headings: 
I. "Effective G-functions" . 

We first calculated the "effective G-functions" defined by 

Qij) ^ <MiM2\0jG{mg,fi,q^)\B> ( ^ < M^M2\0jGim,,, fi,q^)G{mg„ fi,q^)\B > 
^ <MiM2\0j\B> ' '"''^ <MiM2\0j\B> 

(20) 

where Oj are the different currents. The numerator in eq.(20) is calculated accord- 
ing to eq.(12) where the coefficient Cj(g^) for the current product Oj is replaced by 
function G{mg,fi, q^) but the Lorentz structure of the current product Oi is kept. The 
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nia + nib 



M 



in 



(17) 
(18) 
(19) 



denominator is calculated according to eq.(12) with Ci{q^) = 1. The effective G's are 
independent of j. Thus we drop the superindex henceforth. Gg^gj is relevant to calcu- 
lations involving a quark- ant iquark loop insertion in the gluon propagator. We show 
the results in Table I where we also contrast the effective G-function with the value 
calculated in the conventional approach with two chosen values of q^: = ml /A and 
= ml/2. We notice, first, that Gq are process dependent and, second, the largest dif- 
ference between the standard approach and ours is seen for processes involving charm 
quark in the penguin loop. 

II. CP Asymmetry. 

We calculated the CP asymmetry defined as 

r - r 

= _ (21) 

where T represents the rate for the CP conjugate reaction. This calculation was done 
in three different scenarios: (a) without the inclusion of 0{as) quark-antiquark vacuum 
polarization insertion in the gluon propagator, (b) with the 0{as) vacuum polariza- 
tion loop in the gluon propagator, and (c) without the 0{as) vacuum polarization 
loop and with the imaginary part from the diagonal quark — antiquark states (for 
example, cc contribution for channel D^D^ ) removed, as motivated by CPT 

requirement and advocated in P]. Table 2 summarizes the results of these calcula- 
tions: K^K^ decay is not represented in this Table because it is a pure pen- 
guin process whereas the CPT cancellation occurs between Tree® 0{al) Penguin and 
0{as)Penguin ® O {as) Penguin^, [11], ^ . 



III. Comparison with other model calculations. 
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We compared the results of calculating the decay rates and CP asymmetries in four 
different methods as listed below: 

(i) B-S approach: Here we calculate the decay rates and in our approach but with 
the gluon vacuum ploariztion loops and the specelike penguin graphs, Fig. 2(c), 
omitted. 

(ii) As the B-S approach also allows us to calculate the transition form factor -Fo(q'^) 
(defined by Bauer, Stech and Wirbel (BSW) |13[)) using the wavefunctions (14) and 
(15) in a covariant formalism, we attempted a conventional calculation (in the 
factorization approximation ) with the calculated B-S form factor FQ{q^) and Gq{q^) 
evaluated at = ml/4 and 1^1/2. The B-S formfactor FQ{q^) can be written down in 
a form similar to eq. (11-13). 

(iii) A hybrid calculation: we attempted a conventional calculation ( in the 
factorization approximation ) with BSW|]TB[ form factors and the effective G 



function, Gq{q^) of (20) rather than the G-function itself. 

(iv) A conventional calculation (in the factorization approximation) with BSW[0 
form factors and Gq{q^) evaluated at = and ml/ 2. 

The results are shown in Table 3. 

We need to point out the following features of Table 3. First, we found, as was also 
the case in [Q, that the absolute branching ratios came out too low by more than two 
orders of magnitude. Further, the shortfall was decay mode dependent. The branching 
ratios appearing in Table 3 under the columns marked Method 1 and 2 were normalized 
according to the following prescription: 
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Experimental values of B{B^ D l^vi) and B{B^ tt l^vi) are from [14| and 



||15|| , respectively. Despite the need to thus normalize the branching ratios, the CP 
asymmetry ought to be reliable. Second, in our B-S approach, q is integrated over, but 
in methods where q has to be chosen (Methods 2 and 4), acp depends sensitively on 
the choice of q^. This is due to the fact that the imaginary part of the cc penguin loop 
depends strongly on the choice of in the range 

/4 < g2 < ml/2. 

5. Discussion 

We begin by listing the parameter set we have used in our calculations. 

mu = md = 5MeV, m, = 160MeV, = 1.35GeV, rrib = b.OGeV. 
fB = fD = 2QQMeV, /k and U from[lA] 
as = 0.7lGeV, an = OAlGeV, ax = 0.33GeV", = 0.33GeV 



GKM angles : A = 0.90, A = 0.22, p = -0.12, r] = 0.34p| (24) 



In calculating errors in the CP asymmetry acp , we allowed a ±10% error in the above 
listed Ob, C(D and ok- To put this in perspective in terms of the variation of parame- 
ters in (17)-(19), allowing mb to vary in the range (4.8-5.2) GeV and Jb in the range 
(180-220)MeV, ob changes by ±7% around its central value; similarly, allowing rric 
to vary in the range (1.2-1.6) GeV and fo in the range (180-220)MeV, an changes 
by ±12% around its central value; and allowing vary in the range (150-200) MeV 
changes ax by ±7%. 



We have done a calculation of CP asymmetries in a selected few hadronic two-body 
B decays in a Bethe-Salpeter formalism and contrasted the results with those of the 
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conventional factorized approach. We have also investigated the effect of including the 
gluon quark-antiquark vacuum polarization loops. We find that Ocp depends strongly 
on the method adopted to calculate it. In particular, for decays into charmless mesons, 
the conventional factorization approach results in an asymmetry which depends very 
strongly on the choice of (see Methods 2 and 4 of Table 3). This is due to the fact 
that cc loops generate an imaginary part which is sensitive to the choice of in the 
range < q^ < fnl/2. Our method, on the contrary, smears over q^ through a 

convolution of meson wavef unctions. 



A shortcoming of our approach (which we share with [^) is that the branching ratios 
are too low. A single overall normalization, independent of the decay mode, does not 
correct the situation. Thus, we normalized B[B^ D^D^) to the measured branching 
ratio for semileptonic 5° D-l+vi, and B{B- R-n^) and B{B- R-R^) to 
the semileptonic B^ — > tt^I^ui branching ratio. However, a^p ought to be reliable. In 
technical terms the low value of the decay rates is related to the time-oscillators which 
further dampen the contribution of the space-oscillators in our covariant wavefunctions. 
Perhaps the use of more sophisticated bound state wavefunctions IT^ would alleviate 
this difficulty. 
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Figure Captions 

Fig. 1(a): 0{as) penguin graph for b — > sqq transition. 

Fig. 1(b): Penguin graph for b — > sqq transition with 0{as) gluon quark-antiquark 
vacuum polarization loop. 

Fig. 2 (a): Tree graph arising from operators Oi and O2 leading to the amplitude in 
eq.(ll). 

Fig. 2(b): Graph arising from operators O3 — Oiq leading to the amplitude in eq.(12). 
Fig.2(c): Graph arising from operators O3 — Oio leading to the amplitude in eq.(13). 
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Table 1: The effective G functions of Penguin diagrams for some processes. See Eq.(20 
) for definitions. Last two rows list the values of Gq ai, — ml/A and — ml/2. 



Channel 


Gu 


Gs 


Gc 


Guu 


Gus 


Guc 


Gcs 


Gcc 


B- D-D° 


-6.7+1.8i 


-2.1+1.81 


0.24+0.321 


42-251 


11-161 


-4.1-1.41 


-1.9+0.101 


-0.095+0.571 


B- K'tt" 


-7.1+1.4i 


-2.5+1.41 


0.19+0.671 


49-211 


16-141 


-2.7-4.61 


-1.9-1.51 


-0.81+0.411 


B- K'K^ 


-7.1+1.4i 


-2.5+1.41 


0.19+0.681 


49-211 


16-141 


-2.8-4.71 


-1.9-1.51 


-0.82+0.411 


= ml/A 
= ml/2 


G„ = -7.2 + 2.1i Gs = -2.f> + 2.li = 0.87 

G„ = -7.6 + 2.1i Gs = -3.0 + 2.1i Gc = 0.65 + 1.7i 



Table 2: CP asymmetries acp{%) under different scenarios: column 2, without gluon 
vacuum ploarization loops; column 3, with quark- antiquark gluon vacuum polarization 
loops; column 4, without gluon vacuum polarization loops and imaginary parts arising 
from diagonal channels removed. Figures in parentheses correspond to a calculation 
where spacelike penguins (Fig. 2(c)) were ignored. 



Channel 


gluon vacuum 


gluon vacuum 


Imaginary part from 




ploarization excluded 


polarization Included 


diagonal channels removed 


B~ D-D° 


2.7 (2.9) 


2.1 (2.2) 


2.7 (2.9) 


B- K'TT^ 


4.0 (2.9) 


3.1 (2.3) 


3.4 (2.4) 
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Table 3: The calculated CP asymmetries and branching ratios in different methods 
Method 1: 0{as) calculation in our B-S approach; Method 2: Factorization approxi- 

2 2 

mation with B-S form factors. First entry uses (f = ^ and second = Method 
3: Factorization approximation with BSW from factors and effective Gq] Method 4: 

2 

Factorization approximation with BSW form factors. First entry uses = and 
second — 



Channel 


Method 1 
Br and acp{%) 


Method 2 
Br and acp{%) 


Method 3 
Br and acp(%) 


Method 4 
Br and acp{%) 


B- D'D^ 


5.3xl0-^ 2.9 ±0.1 


5.2 X 10"^, 3.3 

5.3 X 10-^, 3.6 


4.7 X 10"^, 2.8 


4.6 X 10"^, 3.2 

4.7 X 10-^, 3.5 


B- K'TT^ 


7.0xl0-^ 2.9 ±0.2 


7.8 X 10-^ 0.68 
8.0 X 10-^ 6.2 


1.1 X 10-^ 4.6 


1.2 X 10-^ 0.91 
1.3 X 10-^ 9.4 


B- K'K^ 


4.0 X IQ-^ -7.1 ±0.6 


4.3 X 10-6, -19 
5.1 X 10-6, -5.3 


1.0 X 10-6, -7.0 


1.1 X 10-6, -19 

1.2 X 10-6, -5.3 
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